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TIME-DEPENDENT  PROPAGATION 
OF  HIGH-ENERGY  LASER  BEAMS 
THROUGH  THE  ATMOSPHERE:  III 


ABSTRACT 

An  accurate  computational  model  of  combined  strong  t  ^-blooming  and  steady-state 
multiple-pulse  blooming  is  described.  A  single-pulse  calculation  and  a  series  of  calculations 
for  an  overlap  number  of  1.5  are  presented.  A  new  method  is  derived  for  calculating  a 
corrective  phase  to  compensate  for  thermal  blooming.  This  method,  which  is  based  on  an 
equivalent  thin  lens  model,  provides  a  fast  and  accurate  estimate  of  the  compensating  phase 
when  most  of  the  phase  distortion  occurs  near  the  laser.  It  is  capable  of  computing  the 
steady-state  compensation  for  thermal  blooming  of  either  cw  or  multiple-pulse  laser  beams 
in  noncoplanar  scenarios.  Several  examples  of  phase  compensation  calculations  are  pre¬ 
sented  including  compensation  for  a  quasi-stagnation  zone  (an  almost  zero  wind-speed 
minimum)  near  the  laser.  Finally,  the  effects  of  realistic  intensity  and  phase  profiles,  stagna¬ 
tion  zones,  and  smoke  clouds  on  laser  beam  propagation  in  Army  scenarios  are  examined. 
In  particular,  simulation  of  the  ABEL  laser  beam  with  realistic  intensity  and  phase  profiles 
shows  that  its  phase  distortion  is  far  more  detrimental  than  its  nonideal  intensity  profile. 
The  stagnation  zone  calculations  suggest  that  the  combination  of  substantial  slewing  rates 
with  a  small  wind  speed  at  the  laser  will  make  stagnation  zones  inconsequential  in  many 
Army  scenarios. 


1.  INTRODUCTION 

This  is  the  third  in  a  series  of  annual  reports  describing  the  development  of  computational  models  of 
high-energy  laser  beam  propagation  through  the  atmosphere.  In  addition  to  various  coplanar  ^  and 
noncoplanar^  steady-state  algorithms,  the  first  two  reports  describe  digital  computer  algorithms  for 
calculating  the  time-dependence  of  thermal  blooming  in  the  isobaric  approximation. 

An  updated  list  of  the  features  of  the  computer  models— collectively  known  as  the  FOUR-D  code— is 
found  in  Table  1.  Four  features  have  been  added  since  Ref.  2: 

•  Thermal  conductivity  for  use  in  modeling  laboratory  scale  experiments; 

•  A  new  multiple  time-step  /  ^-blooming  model; 

•  A  new  predictive  phase-compensation  method  (the  Bradley  and  Herrmann  model  ’  can  also  be  im¬ 
plemented  for  comparison  purposes);  and 

•  A  simplified  method  for  estimating  the  long-time  average  beam  spread  due  to  turbulence. 


Table  1.  Basic  outline  of  current  FOUR-D  propagation  code. 

Variables  x,y,z,t, 

where  x,  y  are  transverse  coordinates  and  z  is  axial  displacement. 

Form  of  propagation  equation  Scalar  wave  equation  in  parabolic  approximation 

^  =  V?  S  +  -  l)a  • 

OZ  L 
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Table  1.  (continued) 


Method  of  solving  propagation  equation 


Hydrodynamics  for  steady-state 
cw  problems 


Transonic  slewing 


Treatment  of  stagnation-zone 
problems  for  cw  beams 


Nonsteady  treatment  of  multipulse 
density  changes 


Metiiod  of  calculating  density  change 
for  individual  pulse  in  train 


Symmetrized  split  operator,  finite  Fourier  series,  fast  Fourier  transform 
(FFT)  algorithm 


=  exp  (-  ^  V?)  «P  (-  x) 


exp 


n 


X  =  *2  (m2  _  1)  . 

Uses  exact  solution  to  linear  hydrodynamic  equations.  Fourier  method 
for  M  <  1.  Characteristic  method  for  M  >  1.  Solves 

bp  bp 

^  ’’y  Po\'  ^1  =  ’ 

^o[\  ^  Vi  =  »  - 

fx  (pi  •  • 

Steady-state  calculation  valid  for  all  Mach  numbers  except  M  - 
Code  can  be  used  arbitrarily  close  to  M  =  1. 

Time-dependent  isobaric  approximation.  Transient  succession  of 
steady-state  density  changes;  i.e.,  solves 


3Pl  3P,  ,  , 

IT  ^  ■'x  3  -  V?  Pi  =  (T  -  1)  acf/. 

X 

Changes  in  density  from  previous  pulses  in  train  are  calculated  with 
isobaric  approximation  using 

9Pl  9Pi  - 

=  -  (7  -  l)ac;2  ^  5(,  _ 

n 

where  T/„(x,y)  is  wth  pulse  fluence.  Density  changes  resulting  from  the 
same  pulse  are  calculated  using  acoustic  equations  and  triangular 
pulse  shape. 

1)  Triangular  pulse  model  for  moderate  -blooming: 

Takes  two-dimensional  Fourier  transform  of 


where  I  is  Fourier  transform  of  intensity,  and  t  is  the  time  duration 
of  each  pulse.  Source  aperture  ^ould  be  softened  when  using  this 
code  provision. 
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Table  1.  (continued) 


Treatment  of  steady-state  multipulse 
blooming 


Treatment  of  turbulence 


Lens  transformation  and  treatment 
of  lens  optics 


Treatment  of  nondiffraction- 
limited  beams 


2)  Multiple  time-point  model  for  strong  t  -blooming: 

Solves  die  full  set  of  time-dependent  linearized  hydrodynamic  equations 
widiin  the  pulse.  Density  is  the  sum  of  this  solution  and  steady-state 
multipulse  contribution  determined  by  fluence  of  current  pulse. 


Previous  pulses  in  train  are  assumed  to  be  periodic  replications  of 
current  pulse.  Solves 

5(t  -  r  )  . 

Pulse  self-blooming  is  treated  as  in  the  nonsteady-state  case. 


3/ 


+  V 


bx 


+  V 


y  by 


oal 


1)  Uses  phase-screen  mediod  of  Bradley  and  Brown  widi  Von  Karman 
spectrum.  Phase  screen  determined  by 


dky  exp  [i{k^x  +  k^)] 


where  a  is  a  complex  random  variable  and  is  spectral  density  of 
refractive  index  fluctuations. 

2)  Alternately,  takes  the  convolution  of  target-plane  intensity  widi  a 
Gaussian  distribution  using  Lutomirski  and  Yura’s  intermediate-range 
beam-spread  formula  [See  H.  T.  Yura,  Appl.  Opt.  10,  2771  (1971)] 
The  jitter  is  increased  by  0.0242  z)^*^  where  is  in 

k  in  cm”^,  and  z  in  cm. 


Compensates  for  a  portion  of  lens  phase  front  with  cylindrical 
Talanov  lens  transformation.  Uses  in  spherical  case 


where  Zf  is  focal  length  of  lens,  is  focal  lengdi  compensated  for  by 
Talanov  ^transformation,  and  is  focal  length  of  initial  phase  front. 


1)  Spherical-aberration  phase  determined  by 

^SA  =  22^  (^2  ^  ^2)2 

or 

2)  phase-screen  method  of  Hogge  et  al.  Phase  determined  as  in 
turbulence,  with 


where  is  correlation  length  and  is  phase  variance. 


3 


Table  1.  (continued) 


Adaptive  lens  transformation 


Selection  of  z-step 


Scenario  capability 


Treatment  of  multiline  effects 


Treatment  of  beam  jitter 


Removes  phase 
2 

^  [a/x.  -  +  e.(x.  -  Oc;)j 

/=! 

throu^  lens  transformation  and  deflection  of  beam.  H^e  Xi  -  x,  X2 
=  averages  are  intensity  weighted,  a/  and  are  calculated  to  keep 
the  intensity  centroid  at  mesh  center,  and  intensity  wei^ted  r.m.s. 
values  of  x  and  y  are  constant  with  z  . 

Adaptive  z-step  selection  based  on  limiting  gradients  in  nonlinear 
contribution  to  phase  and  limiting  the  wind-velocity  change.  Constant 
z-step  over  any  portion  of  range  also  possible. 

General  noncoplanar  scenario  geometry  capability  involving  moving 
laser  platform,  moving  target,  and  arbitrary  wind  direction.  In  coplanar 
case,  wind  can  be  function  of  z  and  t. 


Calculates  average  absorption  coefficient  based  on  assumption  of 
identical  field  distributions  for  all  lines 


Takes  convolution  of  intensity  in  target  plane  with  Gaussian 
distribution: 


X  lix  -  ;c',  y  ~  y')  , 

where  o  is  variance  per  axis  introduced  by  jitter. 

Treatment  of  predictive  phase  compensation  1)  Laser  beam  is  propagated  a  distance  Zj,  witii  steady-state  thermal 

blooming,  then  back  to  the  laser  aperture  in  vacuuo.  Difference 
between  new  and  original  phase  is  an  equivalent  thin  lens  whose 
conjugate  is  used  to  correct  for  the  thermal  blooming.  Useful  when 
blooming  phase  change  is  concentrated  near  laser.  Can  be  iterated. 

2)  Bradley  and  Herrmann  scheme  is  also  implemented  [Appl.  Opt.  13, 
331  (1974)].  Useful  for  coplanar  scenarios  in  which  blooming  is 
concentrated  near  laser. 

Code  output  Isointensity,  isodensity,  isophase,  and  spectrum  contours.  Intensity 

averaged  over  contours.  Plots  of  intensity,  phase,  density,  spatial 
spectrum  along  specific  directions,  etc.,  at  specific  times. 

Plots  of  peak  intensity  and  average  intensity  vs  time. 

Numerical  capacity  when  used  with  Spatial  me^,  64  X  64,  30  sampling  times,  no  restriction  on  number 

CDC  7600  and  restricted  to  internal  of  axial  space  increments. 

memory  (large  and  small  core) 

Problem  zoning  features  Number  of  space  increments  in  x  and  y  directions  must  be  equal  and 

expressible  as  a  power  of  2. 


t  ^-Blooming  Model 

In  Section  2  we  describe  a  new  t  ^-blooming  calculation  method.  In  general  three  phenomena  may  require 
a  fully  time-dependent  calculation:  t  ^-blooming,  transonic  slewing,  and  the  occurrence  of  a  null  wind 
velocity.  Of  these  three,  only  strong  t  ^-blooming  is  never  amenable  to  a  steady-state  analysis.  A  transonic 
zone  along  a  slewed  laser  beam  has  only  a  small  effect  on  the  beam,  at  least  for  ranges  under  10  km;  outside 
the  sonic  region,  steady-state  conditions  prevail.  Most  scenarios  do  not  have  a  true  stagnation  point,  but 
rather  only  a  region  of  low,  but  finite,  transverse  wind  velocities  because  the  wind  velocity  rarely  lies  exactly  in 
the  slewing  plane.  Since  the  minimum  velocity  is  often  an  appreciable  fraction  of  1  m/s  for  military  problems, 
most  military  “stagnation  zone”  problems  can  be  calculated  with  a  noncoplanar  steady-state  wind  field.  The 
calculation  of  moderate  or  strong  t  ^-blooming,  on  the  other  hand,  requires  a  time-dependent  treatment. 

Since  typical  multipulse  CO2  lasers  have  pulse  lengths  under  100  ^us,  interpulse  times  of  over  5  ms,  and 
beam  diameters  under  1.5  m,  the  isobaric  approximation  with  a  series  of  impulse  source  terms  adequately 
describes  the  air-density  change  due  to  all  but  the  current  pulse.  Finding  the  density  change  due  to  the  current 
pulse  requires  solving  the  full  set  of  linearized  hydrodynamic  equations.  Within  each  pulse,  the  new  multi¬ 
ple  time-step  t  ^-blooming  method  self-consistently  solves  these  equations  and  the  Fresnel  wave  equation  in 
the  presence  of  an  air-density  contribution  from  a  periodic  replication  of  the  current  pulse  (multipulse  steady- 
state  overlap  blooming).  This  new  method,  which  is  based  on  a  Fourier  transform,  removes  the  restrictions  of 
our  simpler  (but  faster)  trangular-pulse  approximation  and  allows  the  simulation  of  pulses  whose  shape 
changes  drastically  in  time.  It  also  includes  the  effect  of  finite  wind  velocities  within  the  pulse  when  v  >  0. 1 
(to  our  knowledge  no  other  propagation  code  has  included  this  effect);  such  effects  might  be  significant  when 
multipulse  lasers  are  used  against  Mach  0.7  or  faster  targets  in  crossing  target  scenarios. 


Thin  Lens  Model 

The  extent  to  which  adaptive  optics  can  compensate  for  thermal  blooming  is  being  actively  investigated 
with  experiments  and  computer  calculations.  When  the  thermal-blooming  phase  distortion  occurs  near  the 
laser  aperture,  several  simplified  models  adequately  describe  the  required  phase  correction  at  the  laser  aper¬ 
ture.  For  example,  Bradley  and  Herrmann  ^  have  successfully  used  a  multiple  of  the  line  integral  of  the  initial 
beam  shape  along  the  wind  direction  as  the  phase  compensation  for  thermal  blooming  of  cw  laser  beams  at 
moderate  to  high  slewing  numbers  (A^  =  o)R/vq  =  velocity  change/initial  velocity)  in  coplanar  scenarios.  In 
Section  3,  we  describe  a  more  general  model  of  comparable  simplicity  that  computes  the  phase  compensation 
for  steady-state  multipulse  (when  t  ^-blooming  can  be  ignored)  as  well  as  cw  thermal  blooming  in  both 
coplanar  and  noncoplanar  scenarios.  For  the  case  described  by  Wallace  and  Pasciak,  ^  the  latter  method  gives 
results  comparable  with  those  obtained  from  a  complex  ray-trace  simulation  of  an  ideal  phase-conjugate 
adaptive  optics  systems.  All  of  these  methods,  as  well  as  the  work  of  Brown  et  al.,^  seem  to  suggest  that  the 
on-target  intensity  increase  due  to  adaptive  optics  corrections  tends  to  be  modest — a  factor  of  about  2  is 
typical — even  though  in  some  scenarios  the  improvement  can  exceed  a  factor  of  6. 


Recent  Army  Studies 

Sections  4  to  8  describe  a  series  of  thermal  blooming  studies  that  are  relevant  to  the  Army  HEL  program. 
Section  4  compares  the  performance  of  large  and  small  aperture  systems  for  the  Hughes  ABEL  laser.  It  is 
shown  that  increasing  the  aperture  size  by  a  factor  of  2  significantly  improves  the  performance  of  the  ABEL 
laser  system,  especially  at  long  ranges.  Section  5  describes  the  effects  of  stagnation  zones  in  two  typical  Army 
scenarios.  For  these  scenarios,  thermal  blooming  generally  reduces  the  on-target  irradiance  less  than  tur¬ 
bulence  and  jitter.  In  Section  6,  the  performance  levels  of  equal-power,  equal-outer-diameter,  uniformly- 
illuminated  disk  and  annular  laser  beams  in  the  cw  and  high-overlap-number  multipulse  laser  cases  are  com¬ 
pared;  the  disk  aperture  is  shown  to  suffer  less  thermal  blooming  in  agreement  with  the  results  of  other 
groups.  In  Section  7,  thermal  blooming  due  to  a  smoke  screen  located  near  the  laser  is  discussed.  Finally, 
the  results  of  a  study  of  the  relative  importance  of  thermal  blooming  at  two  potential  White  Sands  Missile 
Range  laser  test  sites  are  described  in  Section  8. 
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Time-dependence  of  strong  t  -blooming  for  a  Gaussian  beam,  (a)  Three-dimensional  plot  of  radial  intensity  profiles  vs  time,  (b)  Successive  radial  intensity  profiles. 


2.  MULTIPLE  TIMESTEP  r  ^.BLOOMING  MODEL 


Previous  t  ^-Blooming  Studies 

Both  computational  ^’6  and  experimental  studies  of  strong  t  ^-blooming  in  a  stationary  absorbing  at¬ 
mosphere  indicate  that  initially  bell-shaped  irradiance  profiles  develop  a  deep  on-axis  depression  late  in  the 
pulse.  A  typical  result^  for  a  square-temporal  pulse  and  a  Gaussian  beam  is  shown  in  Fig.  1.  At  the  suc¬ 
cessively  later  times,  the  pulse  shape  flattens  and  broadens,  ultimately  developing  a  deep  central  depression 
with  a  high-intensity  rim.  At  even  later  times,  the  rim  itself  actually  splits  and  broadens.  Figure  2  shows  the 
saturation  of  the  on-axis  fluence  due  to  the  central  depression  in  the  back  of  the  pulse.  When  t  -blooming  is 
severe,  no  steady-state  model  that  ignores  this  intrapulse  time-dependence  can  adequately  simulate  the  ther¬ 
mal  blooming. 

Our  previous  investigations  have  shown  that  a  simple  single-time-step  triangular-pulse  model  ade- 
quately  describes  t  ^-blooming  prior  to  significant  saturation  of  the  on-axis  fluence.  The  computational  speed 
gained  by  using  only  a  single  timestep  makes  this  model  extremely  attractive  when  it  is  applicable.  However,  it 
oversimplifies  the  physics  of  strong  t  ^-blooming  by  assuming  a  fluence  approximately  proportional  to  the  in¬ 
tensity  at  the  center  of  the  pulse.  ■  ,  e-  •  e  u 

The  perturbation  theory  expressions  currently  in  use  are  even  more  drastic  oversimplifications  ot  the 
physics  of  strong  t  ^-blooming  when  overlap-blooming  is  significant.  They  are  based  on  the  assurnption  that 
the  t  ^-blooming  is  driven  by  the  unbloomed  laser  beam  shape,  consequently  they  overestimate  t  -blooming 
by  ignoring  the  increased  laser  beam  spot  sizes  near  the  target  due  to  phase  changes  induced  by  overlap¬ 
blooming  near  the  laser.  Even  if  overlap-blooming  is  negligible,  the  perturbation  expressions  tend  to  un¬ 
derestimate  the  time  required  for  the  on-axis  fluence  to  saturate  by  ignoring  the  self-moderating  negative  feed¬ 
back  of  self-consistent  t  ^-blooming.  Nonetheless,  the  perturbation  expressions  provide  useful  estimates  for 
t  ^-blooming  irradiance  reductions,  if  overlap-blooming  is  negligible. 


Description  of  Model 

Within  a  representative  pulse,  the  FOUR-D  code’s  new  multiple  time-step  t  ^-blooming  routine  self- 
consistently  solves  the  full  set  of  time-dependent  linearized  hydrodynamic  equations  and  the  Fresnel  wave 
equation: 


-  2ik  ^  ^(r,0  +  Vx  ^(^>0 
2 

+  ^  5p(r,0  a(r,0  =  0  ,  (1) 

"o 

^  +  Hr,f)  =  0  ,  (2a) 

^  +  V?  Pi  ('•.0  =  0  ,  (2b) 


with 


^[pj(r,0  -  cl  Pi(r,0]  =  5(r,0  ,  (2c) 


Time  —  ps 

Figure  2.  Saturation  of  the  on-axis  fluence  for  the  strong 
blooming  case  of  Fig.  1.  The  circles  are  the  results  of  the  triangular 
pulse  t  ^-hlooming  model. 
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where  5p  is  the  total  perturbation  in  the  air  density  including  the  effects  of  heating  by  previous  pulses;  is  the 
ambient  air  density;  Vg  is  the  transverse  wind  velocity  relative  to  the  laser  beam;  p  1,  v  j,  and p  \  are,  respectively, 
the  density,  the  local  flow  velocity,  and  pressure  changes  due  to  the  current  pulse;  /  is  the  instantaneous  inten¬ 
sity  of  the  laser  beam;  F  is  the  single-pulse  fluence;  ^  is  the  slowly  varying  complex  electric  field  amplitude;  At 
is  the  time  between  pulses;  and  the  remaining  quantities  have  their  usual  physical  significance.  The  density 
term  that  contains  F  is  the  multipulse  steady-state  contribution  to  the  air  density  in  which  all  pulses  are 
assumed  to  be  identical.  Equations  (2)  are  solved  with  z  as  a  parameter,  since  the  laser  beam  changes  shape 
only  over  distances  that  are  large  in  comparison  with  its  diameter. 

The  split-step  finite  Fourier  method  of  solving  Eq.  (1)  reduces  the  calculation  to  alternate  free- 
propagation  and  phase-screen  steps.  The  latter  uses  the  density  8p  calculated  from  the  instantaneous  intensity 
and  the  single-pulse  fluence  at  the  location  of  the  phase  screen.  Equation  (2a)  can  be  used  to  eliminate  pj  from 
Eq.  (2c);  this  reduces  Eqs.  (2b)  and  (2c)  to  a  coupled  pair  of  equations  in  the  variables  pi  and  x]/.  After  these 
equations  are  Fourier  transformed  in  the  spatial  variables  x  and  y,  they  are  diagonalized  by  changing  the 
variables: 


• 


(4) 


These  variables,  then,  satisfy 


at 


(5) 


Since?  is  a  known  function  of  t,  the  solution  of  Eq.  (5)  for <  t  <  is  given  by: 


exp 


(6) 


where  (7^  =  ±  .  The  value  of  S  can  be  approximated  by 


5(kj^,r)  =  a^(k^)  (t  -  tj  +  ?„(k^)  for  >  t  >  t^. 
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Then,  at  the  next  time  step  is  given  by 


(k^)  =  exp  (-iq^At)  +  -  t%) 

X  ql^  •  2i  Atj  exp  (-|  iq±  At)  -  ia„Atq~^  .  (7) 


After  0+  and  $1  are  obtained  from  Eq.  (6),  pi  can  be  obtained  from  Eqs.  (4)  and  the  density  p  i  is  given  by  the 
solution  of  Eq.  (2c): 


or  more  explicitly  by 

~«+l  =  _  i  l^n  +  0«^jexp  i-i\\QAt)  + 

-  %  \}^  (^0**^1)”^  •  k^At^  exp  j  iv^  •  k^At^  -i  (Vg  •  A?j 

-  \  [2  (l  ''o  •  K  (- 1  '■''0  •  ‘'i 

Equations  (7)  and  (9)  mathematically  describe  the  intrapulse  density  update.  Three  arrays  corresponding 
to  ,  ?'i ,  and  ^ must  be  carried  from  timestep  to  timestep  (the  latter  two  to  provide  values  of  ^i). 

Previous  numerical  solutions  of  strong  t  ^-blooming  have  ignored  the  wind  velocity  within  the  pulse.  If 
the  effective  wind  velocities  near  the  target  exceed  Mach  0.7,  intrapulse  wind  velocity  effects  will  be  significant 
(if  t  ^-blooming  is  important  at  all).  The  reason  is  illustrated  in  Fig.  3  for  a  Mach- 1.5  wind  across  a  Gaussian 
laser  beam.  The  density  profile  is  no  longer  centered  on  the  laser-beam  axis  (located  at  the  +),  but  is  forced 
downwind.  After  a  sound  transit  time,  it  begins  to  develop  the  Mach  stem  that  characterizes  the  steady-state 
density  profile.  Its  shape  does  not  even  vaguely  resemble  the  V^/  profile  that  is  characteristic  of  zero  wind 
velocity. 


Results  of  Calculations 

As  a  test  of  the  FOUR-D  code’s  new  multiple  timestep  t  ^-blooming  model,  the  single-pulse  calculations 
of  Figs.  1  and  2  for  1.8-km  range  were  repeated.  The  initial  laser  beam  was  assumed  to  have  an  infinite  Gaus¬ 
sian  shape  with  a  50-cm  e  intensity  diameter.  The  beam  was  focused  at  2.5  km.  The  temporal  pulse  shape 
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Figure  3,  Evolution  of  heated  air-density  profiles  in  a  Mach-1.5  wind.  The  scale  lines  are  t,  where  is  the  speed  of  sound.  The  heat  source 
is  a  Gaussian  laser  beam. 

was  assumed  to  be  square,  with  a  133-/iS  duration,  and  the  pulse  energy  was  assumed  to  be  62  kJ.  The  laser 
wavelength  was  10.59  /im,  and  the  atmospheric-absorption  coefficient  was  assumed  to  be  0.3  km  ~K  The  beam 
quality  was  2  times  diffraction-limited,  represented  by  wavelength-scaling  throughout.  The  results  are  shown 
in  Fig.  4.  For  approximately  the  first  10  ms,  the  beam  profile  remains  nearly  Gaussian  as  the  heated  air  has  not 
yet  had  time  to  expand  and  create  a  significant  density  change.  By  33.3  ms,  the  density  change  along  the  laser 
beam  has  become  significant  and  the  laser  intensity  profile  is  noticeably  more  flat-topped  than  the  initial 
Gaussian.  By  66.6  ms,  the  intensity  has  a  pronounced  on-axis  depression  and  a  shape  like  a  doughnut.  At 
100  ms,  the  central  depression  has  become  deeper  and  the  doughnut  radius  has  expanded.  By  133.3  ms,  the 
doughnut  rim  itself  has  widened  and  developed  its  own  central  depression.  The  maximum  intensity  vs  time  for 
this  same  calculation  (Fig.  5)  decreases  steeply  until  the  depression  develops,  plateaus  as  the  rim  merely  ex¬ 
pands,  and  then  drops  again  as  the  rim  is  bloomed.  Qualitative  agreement  of  the  present  model  with  the  results 
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Figure  4.  Single-pulse  (< ^-)  blooming  in  the  Four-D  code.  Time-dependence  of  the  radial  intensity  profiles  under  strong  ?  -hlooming  condi- 
tions. 


Time  —  ms 


Figure  5.  Time-dependence  of  the  maximum  intensity  for  the 
calculation  of  Fig.  4.  The  upper  flat  curve  is  the  unbloomed  inten¬ 
sity. 


of  the  previous  detailed  radially-dependent  code  is  obvious;  the  qualitative  agreement  is  satisfactory  also.  The 
FOUR-D  code  customarily  uses  less  computer  time  than  the  radially-dependent  code.  It  is  also  generally  sim¬ 
pler  to  use  because  of  its  adaptive  lens  transformation,  and  it  can  perform  a  much  voider  class  of  thermal- 
blooming  calculations  because  it  is  not  restricted  to  radially  symmetric  problems. 
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Figure  6  illustrates  combined  overlap-  and  t  ^-blooming  when  Nq  =  1.5,  =  319,  and  Nf  =  29.6.  A 

laser  beam  with  2.0-MW  power,  15-Hz  rep.  rate,  and  a  1-m  truncated  Gaussian  profile  was  simulated;  it  was 
focused  at  3  km  and  the  target  was  at  2.5  km.  At  6.7  and  13.3  /as,  only  overlap-blooming  is  significant;  by 
20  /AS,  strong  t  ^-blooming  is  evident.  Strong  t  ^-blooming  distortion  at  and  beyond  20  /as  does  not  seem  to 
alter  the  overlap-blooming  much,  as  can  be  seen  from  the  similarity  of  the  overlap-only  irradiance  and  the  6.7- 
/AS  irradiance.  This  might  form  the  basis  of  a  simple  scaling  law  for  t  ^-blooming,  in  which  the  overlap¬ 
blooming  is  estimated  first  and  t  ^-blooming  is  estimated  using  the  overlap-bloomed  laser  beam  area. 

Near-vacuum  values  of  the  on-target  peak  time-averaged  intensity  were  obtained  in  spite  of  the  large  dis¬ 
tortion  number.  In  the  calculation  without  t  ^-blooming,  this  intensity  was  38%  larger  than  the  1 1.4-kW/cm  ^ 
on-target  peak  time-averaged  intensity  of  the  same  beam  propagated  through  a  vacuum  but  corrected  for  the 
linear  absorption  power  loss;  this  focusing  effect  when  the  overlap  number  is  between  1  and  2  has  been  noted 
before.  ^  With  t  ^-blooming  included,  the  on-target  peak  time-averaged  intensity  is  72%  of  the  vacuum  value. 


Figure  6.  Time-dependence  of  on-target  irradiance  for  strong  -blooming  combined  with  overlap  blooming  at  -  1.5. 
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The  1 3-blooming  also  forces  the  instantaneous  intensity  centroid  of  the  laser  beam  to  successively  farther 
upwind  positions:  3.1,  3.3,  5.8,  and  7.5  cm  at  successive  time  increments  of  6-2/3  ^lS.  Because  the  powth  of  the 
density  change  during  a  pulse  is  due  to  the  propagation  of  a  sound  wave  driven  by  the  laser  heating  of  the  air, 
the  density  change  lags  behind  the  laser  heating.  The  shape  of  the  laser  beam  at  early  times  is  the  dominant 
factor  in  determining  the  density  seen  by  the  laser  beam  at  late  times.  This  lag  probably  explains  the  motion  of 

the  intensity  centroid  during  the  pulse.  .  ,  ,  .  .  .  ■  r?-  n 

The  power-dependence  of  the  combined  overlap-  plus  t  3-blooming  for  this  case  is  shown  in  Fig.  7.  Ihe 
overlap  number,  repetition  rate,  and  initial  peak  intensity  were  held  fixed,  and  only  the  pulse  length  was 
varied  For  this  method  of  power  adjustment,  overlap-blooming  increases  the  diffraction-hmited  intensity 
more  than  t  ^.blooming  reduced  it,  until  the  power  exceeds  860  kW.  Above  860  kW,  t  3-blooming  dominates 
and  the  on-target  intensity  is  saturated.  Note,  however,  that  this  calculation  used  a  particularly  favorable 

overlap  number-one  leading  to  a  weak  focusing  (plus  deflection)  of  the  laser  beam.  ,  .  u 

Varying  the  power  by  changing  the  repetition  rate  while  maintaining  the  peak  intensity  and  using  t  e 
26.7-ms  pulse  length  (the  2-MW  pulse  length  in  our  calculations)  would  change  the  results  in  Fig.  6  in  the 

•  The  curve  without  t  ^-blooming  would  follow  the  linear  diffraction-limited  behavior  up  to  higher 
power  levels,  joining  the  present  curve  at  2  MW.  It  would,  thus,  have  a  larger  slope  at  2  MW. 

•  The  fractional  decrease  in  the  on-target  intensity  due  to  t  -blooming  would  be  at  least  as  large  at  all 
lower  power  levels  as  it  is  at  2  MW  in  our  calculation,  since  the  pulse  energy  and  pulse  duration  upon  which 
1 3-blooming  depends  would  then  be  constant.  The  two  curves  would  remain  separate  down  to  the  lowest 

Becaus^the  on-target  peak  intensity  at  the  1-MW  power  level  exceeded  20  MW/cm  3,  we  conclude  that, 
for  overlap  numbers  of  1.5  or  less  and  pulse  lengths  under  100  /is,  the  propagation  of  1-m-optics  CO2  laser 
beams  out  to  ranges  exceeding  2.5  km  will  be  limited  primarily  by  the  dirty  air  breakdown  threshold,  rather 
than  by  thermal  blooming.  Multipulse  COj  laser  beams  with  multimegawatt  power  levels  will  have  to  be 
focused  beyond  the  target  to  avoid  breakdown  during  the  first  few  microseconds  of  the  pulse. 


Figure  7.  Average  on-target  intensity  vs  input  power.  The  power 
was  varied  by  altering  the  pulse  length  while  keeping  all  other 


3.  AN  EQUIVALENT  THIN  LENS  MODEL 
FOR  THERMAL-BLOOMING  COMPENSATION 


Introduction 

Several  optical  techniques  have  been  investigated  for  reducing  the  distortion  of  high-energy  laser  beams 
propagated  through  the  atmosphere.  These  include  coherent  optical  adaptive  techniques  (COAT),  which  are 
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based  on  receiving  a  phase  front  in  the  transmitter  aperture  from  an  unresolved  spot  (glint)  on  the  target.  Two 
experimental  forms  of  COAT  have  been  investigated.  The  first  of  these  is  the  phase-conjugate  COAT 
System,  wherein  the  phase  from  the  glint  is  measured  by  interferometric  means,  and  the  conjugate  of  this 
phase  is  in  turn  applied  as  a  correction  to  a  phased  array  of  mirrors.  The  system  is  adaptive  in  that  the 
monitoring  of  the  return  wave  and  the  adjustment  of  the  phased  array  are  carried  out  continuously.  A  second 
COAT  system  is  the  multidither  scheme,  which  relies  on  a  dithering  or  modulation  of  the  elements  of  the 
phased  array  to  maximize  the  intensity  of  the  glint  return.  The  performance  of  both  the  phase-conjugate  and 
the  multidither  systems,  including  the  finite  resolution  of  both  the  adaptive  mirror  and  the  phase  sensor  array, 
has  been  simulated  by  Brown. 

It  is  also  possible  to  determine  a  phase  correction  by  a  purely  predictive  method  that  depends  only  on  a 
knowledge  of  the  relevant  atmospheric  parameters  and  does  not  require  measured  information  from  a  target- 
glint  return.  Since  the  two  COAT  systems  are  adaptive,  they  can  in  principle  correct  for  both  thermal  bloom¬ 
ing  and  atmospheric  turbulence.  Predictive  methods,  on  the  other  hand,  are  restricted  to  thermal-blooming 
compensation,  since  the  relevant  detailed  information  on  atmospheric  turbulence  is  never  available. 

A  predictive  phase-correction  method  applicable  to  cw  thermal  blooming  was  first  demonstrated  by 
Bradley  and  Hermann.  ^  This  method  assumes  the  phase  correction  to  be  proportional  to  the  thermally- 
induced  density  change  at  the  transmitter.  The  constant  of  proportionality  is  initially  estimated  by  integrating 
the  reciprocal  of  the  wind  speed  along  the  propagation  path,  but  it  is  finally  adjusted  by  iteration.  This  method 
works  well  when  the  atmospheric  heating  is  concentrated  near  the  aperture,  but  it  works  better  for  cw  beams 
than  for  multipulse  beams.  This  is  apparently  due  to  the  fact  that  the  density  variation  with  propagation  dis¬ 
tance  cannot  be  extrapolated  as  easily  for  multipulse  beams  as  it  can  for  cw  beams. 

A  second  predictive  scheme  applicable  to  multipulse  laser  beams  is  due  to  Wallace  and  Pasciak.  ^  This 
method  determines  the  conjugate-phase  correction  by  following  ray  trajectories  from  a  point  source  in  the 
focal  plane  back  through  the  heated  medium  to  the  transmitter  in  a  geometric  optics  approximation.  The 
method  also  utilizes  the  time  development  of  the  thermal  blooming  to  calculate  the  phase  correction  and  the 
atmospheric  heating  self-consistently.  This  can  be  done  because  the  thermal  blooming  for  a  given  pulse  de¬ 
pends  solely  on  the  atmospheric  heating  by  previous  pulses  and  can  thus  be  considered  a  strictly  linear  effect. 

In  this  report  we  describe  a  new  predictive  technique  that  is  based  on  a  thin  lens  approximation  of 
thermal  blooming.  It  compares  favorably  in  accuracy  with  the  two  previously  mentioned  methods.  It  also  has 
some  practical  advantages,  including  ease  of  application  and  versatility.  It  represents  a  generalization  of  the 
earlier  Bradley-Herrmann  scheme  in  that  it  can  employ  information  on  the  density  distribution  along  the 
entire  propagation  path  if  necessary.  At  the  same  time  it  avoids  some  of  the  complications  that  can  occur  in 
the  modeling  of  phase- conjugate  systems.  It  is  equally  applicable  to  cw  and  multipulse  beams  and  can  be  ex¬ 
pected  to  work  well  whenever  the  effective  atmospheric  lens  is  concentrated  near  the  transmitter,  a  condition 
that  seems  to  be  necessary  for  the  unqualified  success  of  any  COAT  or  predictive  phase-correction  scheme.  In 
some  marginal  cases,  the  method  may  give  useful  corrections  if  only  a  portion  of  the  atmospheric  lens  near  the 
aperture  is  compensated. 

Because  of  its  simplicity — determination  of  a  phase  correction  requires  only  the  computation  of  the 
Fourier  transform  of  a  propagated  field — the  method  is  readily  applicable  with  adaptive  coordinate  transfor¬ 
mations  and  noncoplanar  geometry,  which  are  required  for  many  practical  analyses.  Noncoplanar  geometry  is 
particularly  useful  for  studying  the  effects  of  stagnation  zones,  All  real-world  scenarios  are  in  some 
measure  noncoplanar,  and  noncoplanarity  provides  a  nonvanishing  vertical  transverse  wind  component 
within  the  stagnation  zone  that  allows  a  greatly  simplified  steady-state  analysis.  ^ 

The  remainder  of  Section  3  begins  with  a  review  of  basic  equations.  Subsequently,  a  wave-optics  method 
for  computing  a  return-wave  conjugate-phase  correction  is  derived;  however,  no  further  use  is  made  of  this 
model.  Thereafter,  the  thin  lens  representation  of  thermal  blooming  is  described  and  used  to  derive  a  simpler 
phase-correction  method.  Finally,  the  second  method  is  compared  with  the  methods  of  Refs.  7  and  8  in  sam¬ 
ple  calculations.  Results  for  a  representative  scenario  involving  a  stagnation  zone  are  also  given. 


Basic  Equations 

The  electric  field  amplitude  &  is  assumed  to  satisfy  the  Helmholtz  equation  in  Fresnel  approximation 

2  /fc  =  ,  (10) 

az  ■‘- 
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where  the  intensity  /  is  given  in  terms  of  &  and  the  absorption  coefficient  a  by  /  =  |  &  |  ^  exp(-az),Vx  =  a  ^/ax  2 
+  a  2/3y  2  is  the  transverse  Laplacian,  and  Se  =  (ae/ap)5p  is  the  hydrodynamically-induced  change  in  permit¬ 
tivity.  The  steady-state  density  change  for  cw  beams  satisfies 

^  ^  _  (T  -  .11  ^  .  (11) 

^  djc  y  dy  -2 
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The  corresponding  steady-state  density  change  for  multipulse  beams,  on  the  other  hand,  is  governed  by  the 
equation 
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which  is  based  on  the  steady-state  assumption  that  all  pulses  in  the  train  have  the  same  spatial  intensity  dis¬ 
tribution  for  a  given  value  of  z.  In  the  most  general  case  of  a  noncoplanar  scenario  the  transverse  wind  velocity 
vector  V  =  (v^,v^)  rotates  in  the  transverse  plane,  and  the  velocity  components  and  are  Vy(z)  linear  func¬ 
tions  of  z.  Methods  for  solving  Eqs.  (11)  and  (12)  are  described  in  Ref.  2. 

The  numerical  solution  for  the  field  S  2"  in  the  target  plane  can  be  expressed  formally  in  terms  of  the  field 
in  the  aperture  plane  by  means  of  the  equation 

SJ  =  ,  (13) 


where  the  Hermitian  operator  is  defined  by  * 


(14) 


Application  of  is  most  easily  carried  out  in  terms  of  Fourier  transforms,  but  it  is,  in  any  case, 
equivalent  to  a  sequence  of  free-space  propagation  steps  through  free  space  alternating  with  the  imposition  of 
the  phase  fronts  0  =  -(fe/2)  StjAzj,  where  St  is  an  average  of  ht  over  the  path  Azj,  and  the  free-space  propaga¬ 
tion  is  governed  by 

a|S  =  V?s.  (‘5) 

bz  ^ 


Return-Wave  Conjugate-Phase  Correction  Algorithm 

A  wave-optics  analog  of  the  conjugate-phase  method  of  Ref.  8  can  be  derived  with  the  help  of  Eqs.  (13) 
and  (14).  We  wish  to  determine  a  function  yp{x,y)  such  that 

.  (16) 

where  is  the  operator  that  propagates  a  field  from  the  transmitter  to  the  target  through  free  space,  i.e.. 
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(17) 


where  zj  is  the  distance  to  the  target  plane.  In  other  words  multiplication  of  the  field  &  ^by  e  and  propaga¬ 
tion  of  this  new  field  through  the  heated  medium  should  give  the  same  result  as  propagating  the  field  &  ^  the 
same  distance  through  free  space.  It  is  assumed,  of  course,  that  ^  ^  has  the  optimum  phase  for  propagating  in 
free  space. 

If  one  applies  f^e  Hermitian  conjugate  of  to  both  sides  of  Eq.  (16),  the  following  result  is 

obtained: 


.  (18) 

From  Eq.  (14),  it  is  evident  that  applying  the  operator  equivalent  to  propagating  a  field  backward 

from  the  target  plane  to  the  transmitter  aperture.  Thus,  the  right  side  of  Eq,  (18)  corresponds  to  propagating 
the  uncorrected  field  S  ^  through  free  space  to  the  target  plane  and  then  backward  through  the  heated  medium 
to  the  transmitter  aperture  again.  Since  the  function  \p  is  in  general  complex,  Re(t/')  should  be  interpreted  as 
the  negative  of  the  phase  correction. 

Equation  (18)  represents  a  wave-optics  generalization  of  the  conjugate-phase  correction  method  of 
Ref.  8,  and  as  such  contains  the  effects  of  diffraction.  However,  application  of  Eq.  (18)  does  not  produce  a 
self-consistent  heating  pattern  unless  the  beam  is  iteratively  recycled.  The  convergence  or  nonconvergence  of 
such  iterative  schemes  have  already  been  examined. 

Equation  (18)  is  straightforward  to  implement  with  the  existing  propagation  machinery  in  a  wave-optics 
propagation  code.  The  only  requirement  is  that  the  phase-front  information  from  the  nonlinear  propagation 
calculation  be  saved.  Return-wave  schemes  like  the  one  above  have  one  practical  drawback,  however.  In  cases 
of  severe  thermal  blooming,  it  is  useful  to  remove  rapidly-varying  phase  fronts  from  the  beam  by  means  of  an 
adaptive  transverse  coordinate  transformation.  The  required  phase-front  information  will  thus  be  stored 
on  a  set  of  contracting  coordinate  grids.  The  return-wave  calculation  may  have  significant  energy  on  the 
boundaries  of  these  grids,  causing  the  numerical  calculation  to  abort.  The  algorithm  derived  in  the  following 
section  circumvents  this  possibility  while  preserving  some  of  the  wave-optics  features  of  the  method  just 
discussed. 


Thin  Lens  Approximation  of  Heated  Atmoshperic  Lens 
and  Free-Space  Return-Wave  Algorithm 

Let  us  assume  that  the  effect  of  the  distributed  atmospheric  lens  can  be  represented  by  a  thin  lens  placed 
at  the  transmitter  aperture.  We  may  then  write 

•  09) 

Equation  (19)  is  a  good  assumption  if  the  heating  of  the  atmosphere  by  the  laser  beam  is  concentrated  near  the 
aperture.  In  any  case,  one  can  always  attempt  to  represent  by  a  thin  lens  the  density  changes  along  a  limited 
portion  of  the  propagation  path  near  the  aperture,  using  Eq.  (19).  In  that  case,  zjin  Eqs.  (14)  and  (17)  should 
be  replaced  by  the  value  27^,  where  0  <  27^  <  z^,  and  2/  is  to  be  determined  by  trial. 

Applying  the  Hermitian  conjugate  operator  to  both  sides  of  Eq.  (19)  gives 

g/V-gO  =  (20) 

The  right  side  of  Eq.  (20)  describes  the  propagation  of  the  uncorrected  beam  over  a  distance  zj  through  a 
heated  medium  followed  by  propagation  of  the  resulting  field  backward  through  free  space  to  the  transmitter. 
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Computation  of  ^  thus  requires  a  single  Fourier  transform  and  no  storage  of  phase  fronts.  Since  a  thermally- 
bloomed  beam  is  larger  than  the  corresponding  diffraction-limited  beam,  the  return  beam  is  usually  smaller 
than  the  original  uncorrected  beam  and  easily  fits  inside  the  aperture  grid.  Determination  of  a  best  correction 
requires  varying  the  parameter  zj'.  When  the  optimum  value  of  zy  has  been  determined,  an  improved  correc¬ 
tion  may  be  obtainable  by  iterating  to  find  a  self-consistent  heating  distribution.  However,  experience  has 
shown  that  such  additional  iterations  usually  result  in  only  marginal  improvement. 


Numerical  Examples 

The  first  numerical  example  has  been  treated  in  Ref.  7.  The  pertinent  data  are 


Np  =  ka^lR  =  18.3 
=  ocR  =  0.4 


where  Nf  =  Fresnel  number,  =  absorption  number,  =  slewing  number,  and  N[)  =  distortion  number. 
Here  a  is  the  1  /e-radius  of  a  Gaussian  beam  that  has  been  truncated  at  its  1/e  -radius,  and  the  range  R  -  Zfis 
taken  equal  to  the  initial  radius  of  curvature  of  the  beam  zy.  The  isointensity-contour  plots  for  the  corrected 
and  the  uncorrected  beam  are  displayed  in  Figs.  8a  and  8b.  The  correction  was  based  on  Eq.  (20).  A  phase 
correction  ^<t>  based  on  the  method  of  Ref.  7  was  also  calculated: 


A<p 


-ya^N^ 


ln(l  +  NJ 


|&(z  =  oy.yyp’  dx'  . 


(21) 


Computations  based  on  the  two  phase  correction  methods  are  compared  in  Table  2. 

In  Table  2  the  intensities  have  been  averaged  over  the  isointensity  contour  containing  one-half  the  beam 
power.  The  value  of  7  in  Eq.  (21)  and  the  ratio  zj/R  have  been  optimized  to  yield  the  maximum  value  of 


Figure  8.  Isointensity  contours  for  cw  beam  in  focal  plane.  Con¬ 
tours  decrease  from  peak  value  in  increments  of  0.1.  Problem 
parameters  are  18.3,  =  0.4,  =  8,  =  107 ;  (a)  un¬ 

corrected  beam,  (b)  corrected  beam  using  Eq.  (11)  and  Zy//?  = 
0.5.  Correction  decreases  half-power  area  by  a  factor  of  2. 


Table  2.  Comparison  of  results  using  Eq.  (20) 
and  Eq.  (21). 

Eq.  (20) 

Eq.  (21) 

I  //® 

2.09 

1.78 

7 

- 

1.0 

z'plR 

0.5 

- 
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hv/lAV^  ^here  I^y  is  the  average  intensity  for  the  uncorrected  beam.  The  optimum  value  of  IavI^av  for 
Eq.  (20)  is  about  17%  larger  than  the  corresponding  value  for  Eq.  (21).  Further  improvement  in  the  correction 
resulting  from  Eq.  (21)  can  be  obtained  by  readjusting  the  focal  position  Zf,  The  same  presumably  holds  true 
for  Eq.  (20).  However,  no  attempt  has  been  made  to  apply  a  refocusing  correction  here. 

The  second  illustrative  example  is  taken  from  Ref.  8  and  deals  with  a  multipulse  beam.  The  relevant 
parameters  are 


Nq  =  =  4 

^a,=  0 


Na  =  0.4 
Nj)  =  17.9 


where  is  the  overlap  number.  The  beam  is  assumed  to  be  an  infinite  Gaussian  with  1  fe  ^  radius  equal  to  a. 

The  uncorrected  and  corrected  steady-state  isointensity  profiles  at  the  target  are  shown  in  Figs.  9a  and 
9b.  The  uncorrected  peak  Intensity  is  0.52  times  the  value  for  the  diffraction-limited  beam.  Following  the 
correction  based  on  Eq.  (20),  the  peak  intensity  is  0.95  times  the  diffraction-limited  value,  in  essential  agree¬ 
ment  with  the  result  in  Ref,  8.  For  this  example  the  optimum  value  of  zj/R  was  equal  to  1.  An  additional 
iteration  resulted  in  inconsequential  improvement. 

The  final  example  Involves  a  cw  beam  propagating  through  a  noncoplanar  stagnation  zone.  The  horizon¬ 
tal  wind  speed  as  a  function  of  axial  distance  is  seen  in  Fig.  10a,  and  the  velocity  magnitude  as  a  function  of  z 
is  shown  in  Fig.  10b.  The  stagnation  point  (position  of  vanishing  of  horizontal  component  of  wind  velocity)  is 
located  75  m  from  the  aperture.  The  defining  parameters  for  the  problem  are  the  following: 


=  0.183 

=  86.7 

^Z,(0) 

=  450 

=  851 

_  IK^)I 
IKz£,)I 

=  0.727  ] 

Figure  9.  Isointensity  contours  for  multipulse  beam  in  focal  plane. 
Problem  parameters  are  =  9.08,  =  0.4,  =  0,  A^^  == 

18.6.  Initial  beam  shape  is  Gaussian;  (a)  uncorrected  beam,  (b) 
beam  corrected  using  Eq.  (11)  and  z  j/ R  =  1.  For  the  uncorrected 
beam  the  peak  intensity  is  0.53  times  the  peak  for  a  diffraction- 
limited  beam.  Following  correction  this  rises  to  0.95  times  the  dif¬ 
fraction-limited  peak. 
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The  beam  at  z  =  0  is  a  uniformly-illuminated  circular  aperture,  and  the  slewing  number  refers  to  the  stag¬ 
nation  point  z  -  zj).  The  distortion  numbers  are  calculated  at  both  the  aperture  and  the  stagnation  point, 
where  v  ^  ,  the  vertical  component  of  wind  velocity  in  the  transverse  plane,  is  substituted  for  vq.  The  distortion 
and  slew  numbers  are  given  only  to  characterize  the  problem.  Due  to  the  complexity  of  flow  in  the  stagnation 
zone,  these  numbers  cannot  be  used  to  infer  the  thermal-blooming  characteristics  from  those  of  a  coplanar 
scenario  with  similar  values. 

The  uncorrected  and  corrected  isointensity  contours  in  the  target  plane  are  shown  in  Figs.  1  la  and  1  lb, 
and  the  isophase  contours  for  the  corrective  phase  that  must  be  applied  on  the  aperture  mirror  are  shown  in 
Fig.  12. 

The  areas  of  the  uncorrected  beam,  the  corrected  beam,  and  a  diffraction-limited  beam  in  the  target  plane 
are  compared  in  Table  3.  The  area  used  is  that  of  the  isointensity  contour  containing  half  of  the  beam  power. 
Although  the  thermal  blooming  is  clearly  severe.  Table  3  indicates  that  a  phase  correction  can  give  an  im¬ 
provement  by  a  factor  of  about  2  in  the  average  intensity.  This  example  suggests  that  phase  compensation  for 
thermal  blooming  can  be  of  some  use  in  stagnation-zone  situations,  if  the  stagnation  zone  is  situated  near  the 
aperture. 


Figure  10.  Transverse  wind  velocity  as  a  function  of  axial  distance  for  noncoplanar  scenario  involving  a  stagnation  zone, 
(a)  Horizontal  component  vanishes  75  m  from  aperture,  (b)  Transverse  velocity  magnitude  does  not  vanish  at  stagnation  point  due  to 
0.727-m/s  vertical  component. 
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Figure  11.  Isointensity  contours  at  focus  for  stagnation  zone 
scenario:  (a)  uncorrected,  (b)  corrected.  Correction  increases 
average  intensity  by  factor  of  2. 


Figure  12.  Isophase  contours  for  corrected  beam  for  stagnation 
zone  scenario.  Adjacent  contours  represent  phase>change  incre* 
ments  of  7r/4. 


Table  3.  Comparison  of  thermally  bloomed  beam  areas  for 


stagnation  zone 

scenario. 

Unconected 

Corrected 

Diffraction 

limited 

2 

Area  (cm  ) 

12.0 

6.06 

1.71 

z’j-lR 

- 

0.3 

- 

Summary  and  Conclusion 

We  have  described  a  simple  method  for  determining  phase  corrections  to  compensate  for  thermal  bloom¬ 
ing.  The  scheme  is  comparable  in  effectiveness  to  other  predictive  methods,  and  it  is  simpler  to  implement 
than  the  phase-conjugate  return-wave  method  because  representative  density  profiles  need  not  be  saved.  It  can 
be  equally  well  applied  to  cw  and  multipulse  beams  and  it  can  be  used  with  noncoplanar  scenarios  involving 
stagnation  zones. 


4.  MODELING  LASER  PERFORMANCE  WITH 
REALISTIC  INTENSITY  AND 
PHASE  PROFILES:  THE  ABEL  LASER 

The  FOUR-D  code  has  been  modified  to  accept  as  input  an  arbitrary  initial  complex  electric  field.  As  a 
result,  it  can  now  compute  the  performance  of  real  laser  systems  once  their  output  intensity  and  phase  profiles 
are  specified.  The  calculation  of  the  performance  of  the  ABEL  CO2  laser  with  two  aperture  sizes,  small  and 
large,  is  a  good  illustration  of  this  capability.  The  intensity  and  phase  data,  as  well  as  the  propagation 
parameters,  were  supplied  by  Joel  Shafran.  The  propagation  parameters  correspond  to  two  scenarios  of  in¬ 
terest  to  the  Army:  crossing  target  (fly-by)  and  incoming  target  (approach).  Absorption  without  saturation  or 
kinetic  cooling  was  assumed.  Low  frequency  jitter  of  2.39  ^urad/axis  was  modeled.  Turbulence,  on  the  other 
hand,  was  ignored  because  it  is  certainly  negligible,  if  is  less  than  10  "^^m  due  to  the  large  diffraction 
and  thermal-blooming  effects  in  these  calculations. 

The  initial  intensity  and  phase  profiles  of  the  small  (a)  and  large  (b)  aperture  ABEL  laser  systems  are  dis¬ 
played  in  Figs.  13  and  14.  The  sudden  27r-decrease  (one- wave)  on  passing  into  the  oblong  region  on  the  left 
side  of  both  phase  profiles  is  caused  by  a  coding  error  in  the  Hughes  computer  program  that  generated  our  in¬ 
put  data.  Because  the  phase  is  modeled  modulo  27r  in  a  Fourier  code,  this  sudden  phase  jump  has  been 
removed  in  our  calculations,  A  quadratic  term  that  corresponds  to  focusing  on  the  target  was  added  to  the 
phase  in  Fig.  14. 

Our  results  for  the  ABEL  systems  are  summarized  in  Table  4.  From  these  results,  it  is  clear  that  the  large 
aperture  outperforms  the  small  aperture  in  all  these  cases;  for  the  lower  slewing  numbers  in  the  approach 
scenario,  the  improvement  is  substantial.  Such  improvement  is  expected  whenever  the  thermal  blooming  is 
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(a) 


(b) 


Figure  13.  Isointensity  contours  at  the  laser  aperture:  (a)  small  diameter  system,  (b)  large  diameter  system. 


(a)  (b) 


Figure  14.  Isophase  contours  at  the  laser  aperture:  (a)  small  diameter  system,  (b)  large  diameter  system.  The  contour  interval  is  a-/4  (X/8). 
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Table  4.  Relative  performance  of  the  small  and  large  aperture  ABEL  laser  systems. 


Aperture 

and 

scenario^ 

Fresnel 

No. 

(Np) 

Absorption 

No. 

(Na) 

Distortion 

No. 

Slewing 

No. 

Overlap 

No. 

Normalized  on- target  areas^ 

Without  jitter 

With  jitter 

SmaU,  XT 

24.2 

0.45 

825 

301 

87.5 

12.4 

14.3 

Large,  XT 

111.1 

0.45 

415 

301 

188.0 

9.8 

11.2 

SmaU,  IT 

72.6 

0.15 

275 

91 

87.5 

1.5 

1.7 

Large,  IT 

333.4 

0.15 

138 

91 

188.0 

-  1.0 

1.2 

SmaU,  IT' 

24.2 

0.45 

825 

31 

87.5 

107.0 

107.0 

Small,  IT 

24.2 

0.45 

412 

16 

43.8 

43.0 

48.0 

Large,  IT 

111.1 

0.45 

415 

31 

188.0 

14.7 

17.3 

^XT  =  crossing 

target  flying  into  wind.  IT  = 

incoming  target, 

0.1  X 

(XT  range)  offset. 

laser-slewed  into  wind. 

'^Smallest  area 

that  contains 

half  the  power. 

^Marginal  calculation,  probably  overestimates  thermal  blooming. 


large  and  the  overlap  number  for  the  smaller  aperture  is  larger  than  2  (in  these  calculations,  the  overlap  num¬ 
ber  at  the  laser  aperture  is  87.5  for  the  small  aperture).  Figures  15  to  17  are  the  on-target  irradiances.  In  the 
crossing-target  scenario  (Fig  15)  the  on-target  irradiance  has  a  well-defined  central  peak  with  a  broad  wing 
that  is  more  elongated  along  the  direction  perpendicular  to  the  wind  than  parallel  to  it,  and  the  small  aperture 
case  has  considerably  more  substructure  in  the  form  of  secondary  peaks  than  does  the  large  aperture.  The  on- 
target  irradiance  profiles  at  the  shortest  range  in  the  incoming-target  scenario  (Fig.  16)  are  almost  identical  to 
the  corresponding  crossing-target  profiles  (Fig.  15),  except  for  a  scale  factor;  the  intensities  are  about  10  times 
larger  and  the  spot  diameter  is  about  0.4  as  large.  At  the  longest  range  in  the  incoming-target  scenario,  a  long 
ridge  structure  is  the  most  prominent  feature  of  the  irradiance  profiles.  This  elongation  perpendicular  to  the 
direction  of  the  wind  will  probably  become  even  more  pronounced  at  the  longer  ranges  due  to  the  transverse 
defocusing  action  of  the  nearly  stagnant  air  near  the  laser  aperture,  where  the  laser  beam  is  roughly  circularly 
symmetric. 


Figure  15.  On-target  isointensity  contours  for  the  crossing  target  scenario:  (a)  small  diameter  system,  (b)  large  diameter  system.  The 
reference  lines  represent  equal  lengths. 


Figure  16,  On-target  isointensity  contours  for  the  incoming  target  scenario  at  close  range  {Np  -  72.6) :  (a)  small  diameter  system,  (b)  large 
diameter  system.  Both  reference  lines  correspond  to  one-tenth  the  lengths  represented  by  the  reference  lines  in  Fig.  15. 


(a) 


(b) 


Figure  1 7.  On-target  isointensity  contours  for  the  incoming  target  scenario  at  the  longest  range  (Np  -  24.2):  (a)  small  diameter  system,  (b) 
large  diameter  system.  Both  reference  lines  represent  the  same  lengths  as  those  in  Fig.  15. 


All  but  one  of  the  calculations  are  high-quality  computations.  Only  the  largest-distortion-number 
crossing -target  computation  is  less  reliable  due  to  spectral  aliasing.  To  establish  a  lower  bound  to  the  bloom¬ 
ing  in  this  case,  we  repeated  the  calculation  with  twice  the  wind  speed.  As  shown  in  Table  4,  its  area  without 
jitter  is  59%  smaller  than  that  of  lower  wind  speed. 

Comparison  of  the  performance  of  the  actual  to  that  of  more  idealized  laser  beams  shows  that  phase 
imperfections  on  the  ABEL  beams  substantially  reduce  the  on-target  intensities.  In  the  crossing  target 
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scenario  where  the  thermal  blooming  is  smallest,  replacing  the  actual  ABEL  laser  phase  profiles  with  the 
ideal  uniform  phase  increases  the  average  on-target  intensities  for  the  small  and  large  apertures  by  factors  of 
2.2  and  21,  respectively;  in  the  incoming  target  scenario  at  the  longest  range  (where  the  thermal  blooming  is 
greatest),  the  increase  in  on-target  intensity  is  a  more  modest  factor  of  1.8. 

The  nonuniform  intensity  profile  of  the  ABEL  laser  beams  seems  to  have  a  relatively  minor  effect  on  their 
performance.  With  the  ideal  constant-phase  profile,  replacing  the  large  aperture  ABEL  laser  intensity  profile 
with  an  equal  power,  equal  inner-  and  outer-diameter  uniform-annulus  intensity  profile  increased  the  on- 
target  intensity  by  only  21%  in  the  crossing-target  calculation,  where  the  blooming  is  the  weakest.  Further¬ 
more,  in  this  case,  the  uniform  annulus  on-target  intensity  is  only  15%  below  the  diffraction-limited  intensity 
for  a  uniformly-illuminated  circular  aperture  of  the  same  outer  diameter. 


5.  STAGNATION  ZONES  IN  TWO  ARMY  SCENARIOS 


Introduction 

When  the  propagation  path  from  a  high-energy  laser  (HEL)  to  its  target  contains  a  region  of  negligible 
transverse  flow  velocity,  thermal  blooming  can  seriously  degrade  the  on-target  intensity.  For  this  reason, 
stagnation  zones  l»2,5,20  Qf  considerable  interest  within  the  HEL  community.  In  the  relatively  rare, 

but  widely  assumed,  coplanar  scenario  (where  the  laser  beam,  the  target  trajectory,  and  the  transverse  wind 
velocity  lie  in  a  fixed  plane  when  observed  in  the  laser’s  rest  frame),  the  transverse  wind  velocity  vanishes  at  a 
stagnation  point  and  thermal  blooming  attains  a  steady-state  condition  only  when  natural  convection  or 
longitudinal  air  flow  is  considered.  Both  of  these  processes  take  1  to  10  s  to  become  effective  for  typical 
scenarios.  In  the  time  before  these  processes  become  important,  the  phase  shift  due  to  the  region  near  the 
stagnation  point  in  which  the  transverse  flow  has  not  yet  swept  the  heated  air  across  a  beam  diameter  the  so- 
called  stagnation  zone — grows  only  logarithmically  in  time,  so  that  a  quasi-steady-state  develops  even  before 
natural  convection  and  longitudinal  air  flow  establish  a  true  steady-state.  Neverthless,  thermal  blooming  in 
the  presence  of  a  stagnation  zone  between  the  laser  and  its  target  often  severely  limits  the  on-target  intensity. 

In  the  more  general  noncoplanar  scenario,  the  transverse  flow  speed  is  small  (but  nonzero)  at  the  stagna¬ 
tion  point,  which  is  now  defined  as  the  point  on  the  infinitely  extended  laser  beam  where  the  transverse  flow 
speed  is  a  minimum.  Since  the  minimum  transverse  flow  speed  is  often  at  least  an  appreciable  fraction  of 


Figure  18.  The  scenario  for  stagnation  zone  calculations.  The 
horizontal  plane  is  shaded  and  contains  the  wind-velocity  vector. 


1  m/s,  a  steady-state  exists  on  a  timescale  of  the  order  of  <  1  s,  even  when  natural  convection  and  longitudinal 
air  flow  are  ignored.  This  greatly  simplifies  the  calculation  necessary  for  a  realistic  appraisal  of  the  effect  of  a 
stagnation  zone  that  lies  between  the  laser  and  its  target. 

To  evaluate  the  potential  degradation  of  an  HEL  beam  due  to  a  stagnation  zone  in  Army  applications,  we 
have  calculated  the  thermal  blooming  in  a  set  of  scenarios  developed  by  Gebhardt.  The  basic  assumptions 
of  Gebhardt’s  scenarios  are  illustrated  in  Fig.  18.  The  offset,  altitude,  and  speed  of  the  target,  as  well  as  the 
wind  velocity,  are  constant,  and  the  laser  is  stationary;  the  target  altitude  is  half  the  shortest  offset.  The  wind 
velocity  is  2.3  m/s  and  the  absorption  numbers  at  the  longest  range  are  0.19  for  the  DF  laser  and  0.57  for  the 
CO2  laser;  these  approximate  the  summer  (worst  season)  meterological  conditions  at  the  White  Sands  Missile 
Range  HELSTF  site.  The  slewing  number  at  the  longest  range  and  for  the  largest  offset  is  103,  where  the  slew¬ 
ing  number  is  defined  by  the  expression 

[(IWl  -  l''minl)/l''(z  =  0)|]. 

The  laser  beam  was  a  uniformly-illuminated  circular  aperture  and  was  focused  on  the  target.  At  the 
longest  range  and  for  the  largest  offset  the  Fresnel,  overlap,  and  distortion  numbers  for  the  laser  beams  were 
Nf  =  29.9,  No  =  38.5,  and  No  =  650  for  the  CO2  laser;  and  Af  =  83.4  and  Np  =  605  for  the  DF  laser. 


Results 

Thermal  Blooming 

The  main  results  of  our  calculations  are  summarized  in  Table  5.  The  stagnation  point  is  defined  as  the  ax¬ 
ial  position  where  the  transverse  wind  speed  is  a  minimum.  The  stagnation  region  is  defined  to  be  the  propaga¬ 
tion  path  segment  over  which  the  transverse  wind  speed  is  less  than  twice  the  smallest  value  above.  As  il¬ 
lustrated  in  Fig.  19,  low  transverse  wind  speeds  are  confined  to  a  short  region  near  the  laser  in  all  these 
calculations.  In  some  of  these  cases  the  exact  stagnation  point  is  located  slightly  behind  the  laser. 

We  can  draw  the  following  conclusions  from  Table  5: 

•  The  DF  laser  outperforms  the  CO2  laser  in  all  these  cases.  Although  its  strehl  ratio  is  always  lower 
than  the  CO2  laser’s,  the  area  of  its  vacuum  focal  spot  is  ~7.8  times  smaller  than  the  CO2  laser’s,  which  more 
than  compensates  for  its  lower  strehl  ratio.  Additionally,  since  the  CO2  laser  has  less  power  at  a  given  range 
due  to  the  high  absorption  coefficient  at  this  wavelength,  its  performance  is  even  worse  than  the  strehl  ratio 
indicates. 


Table  5.  Stagnation  zone  parameters  and  on-target  intensities. 


Stagnation  region 


Trajectory 

offset 

Target 

location 

Min. 

wind 

speed, 

m/s 

Lengtii 

X  103 

Location® 

X  103 

Av  intensity, 

Strehl  ratio 

CO^  laser  DF  laser 

Nomialized  av 
intensity 

CO^  laser  DF  laser 

0.5 

0 

1.78 

3.15 

0 

0.98 

0.98 

27. 

229. 

0.5 

0.767 

3.38 

0 

1.00 

0.91 

15 

116. 

1. 

0.727 

12.3 

3.39 

0.95 

0.77 

5.6 

41. 

2. 

0.727 

52.0 

22.5 

0.73 

0.35 

1.1 

5.4 

3. 

0.727 

121.0 

57.1 

0.34 

0.151 

0.2 

-  1.0 

1.0’’ 

0. 

1.67 

5.30 

0 

1.01 

0.90 

7.4 

60. 

0.5 

0.869 

3.54 

0 

0.99 

0.90 

5.8 

48. 

1. 

0.400 

3.81 

0 

0.95 

0.79 

3.4 

26. 

2. 

0.394 

20.0 

9.78 

0.84 

0.50 

1.06 

6.4 

3. 

0.394 

40.4 

28.5 

0.51 

0.26 

0.27 

1.6 

®For  those  cases  indicated  by  a  zero,  the  actual  stagnation  point  lies  a  short  distance  behind  the  laser. 
**A11  lengths  are  given  relative  to  this  offset. 
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Figure  19.  Transverse  wind  speed  for  the  shortest  offset  and  the 
longest  range. 


Distance  from  laser  —  fraction  of  range 


•  Even  at  ranges  beyond  the  longest  investigated  here,  the  CO2  laser  would  have  to  be  defocused  to 
avoid  peak  intensities  above  10  MW/cm 

Two  typical  on-target  intensity  profiles  are  shown  in  Fig.  20.  Only  four  of  the  cases  summarized  in 
Table  5  had  intensity  profiles  similar  to  Fig.  20(a):  the  longest  range  DF  laser  cases  in  each  of  Tables  5(a)  and 
5(b),  the  next  to  the  last  DF  laser  case  in  Table  5(a)  and  the  longest  range  CO2  laser  case  in  Table  5(a).  The 
remaining  cases  had  intensity  profiles  similar  to  Fig.  20(b).  Only  the  longer  range  cases  with  their  slower  slew¬ 
ing  rates  had  enough  thermal  blooming  to  significantly  elongate  the  on-target  irradiance. 

-Blooming 

To  estimate  the  decrease  in  the  on-target  intensity  of  the  CO2  laser  due  to  the  t  ^-blooming,  we  repeated 
the  longest  range  calculation  in  Table  5(a),  this  time  including  triangular  pulse  t  ^-blooming.  *  To  simplify  the 


Figure  20.  Typical  on-target  isointensity  contour  plots  (without  jitter  or  turbulence):  (a)  DF  laser,  shortest  offset,  longest  range;  (b) 
CO  2  laser,  longest  offset,  longest  range. 
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numerical  calculation,  the  initial  beam  shape  was  assumed  to  be  an  infinite  Gaussian  beam  of  the  same  power 
and  with  a  1  /e-diameter  equal  to  that  of  the  circular  aperture.  (This  would  overestimate  the  effect  of  t 
blooming.)  The  on-target  irradiance  contours  in  Fig.  21  demonstrate  that  t  ^-blooming  is  minimal;  the  peak 
irradiance  is  reduced  by  only  5%  by  t  ^-blooming.  With  defocusing  to  keep  the  peak  intensity  below  10 
MW /cm  t  ^-blooming  should  be  negligible  at  all  ranges  in  these  scenarios. 

The  above  results  demonstrate  the  danger  in  estimating  t  ^-blooming  from  the  vacuum  propagation  of  a 
laser  beam  when  overlap  pulse  blooming  is  significant.  With  a  perturbation  analysis  that  starts  from  vacuum 
propagation,  Ulrich  and  Hayes  predict  that  t  ^-blooming  reduces  the  on-axis  intensity  at  the  target  to  zero  in 
a  saturation  time  given  in  Eq.  (39)  of  Ref.  2: 


where  N  is  the  refractivity ,  F/j  is  the  pulse  energy,  Tp  is  the  pulse  duration,  and 

is  the  vacuum  beam  area  within  the  1  /e  intensity  contour  for  focal  length/.  The  8.8-^s  value  of  Tj  for  the  above 
problem  is  inconsistent  with  the  minimal  t  ^-blooming  that  we  calculate.  If  we  replace  with  the  area  of  the 
laser  beam  including  repetitive  pulse  blooming,  is  44  fis,  which  is  consistent  with  our  results  of  minimal  t  ^- 
blooming  for  our  20-ms  laser  pulse.  Thus,  overlap-blooming  must,  in  general,  be  considered  for  when  es¬ 
timating?  ^-blooming. 

Jitter  and  Turbulance 

The  effect  of  10  rad  (std.  dev.:  7. 1  rad/axis)  of  jitter  and  turbulence  characterized  by  Cfj  =  10  m  is 
indicated  in  Table  6.  The  turbulence  contribution  to  the  laser  beam  area  is  calculated  from  the  formula  of 
Lutomirski  and  Yura 
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Table  6. 

Contributions  to 

on-target  laser  beam 

area. 

Areas  (cm^) 

Trajectory 

offset 

Target 

location 

Turbulence 

Diffraction 

and 

blooming 

Range 

Jitter 

DF 

CO2 

DF 

CO2 

0.5 

0 

0.559 

0.982 

0.183 

0.121 

0.476 

0.5 

IJl 

0.468 

0.845 

1.0 

1.146 

4.13 

1.82 

1.21 

2.07 

2.0 

13.6 

12.2 

8.09 

2.38 

8.87 

3.0 

29.3 

41.8 

27.7 

12.0 

41.6 

1.0 

0 

3.34 

1.30 

0.860 

0.228 

1.58 

1.132 

4.03 

1.75 

1.16 

0.283 

1.99 

1.0 

1.436 

6.48 

3.74 

2.48 

0.507 

3.25 

2.0 

2.25 

15.9 

15.8 

10.4 

1.95 

9.07 

3.0 

3.172 

31.6 

47.3 

31.4 

7.54 

29.9 

•A.  rjp 

=  nrl=  0.152 

16/5  ^2/5 

(23) 

At  the  longest  range  investigated  here,  the  ratio  of  the  increase  in  the  laser  beam  area  due  to  10  jurad  of  jitter  to 
the  area  that  is  due  to  diffraction  and  thermal  blooming  alone  is  2  to  4  for  the  DF  laser  and  approximately  1 
for  the  CO2  laser.  At  this  range,  10  ^rad  of  jitter  and  10  m  of  turbulence  increase  the  laser  beam  area 
by  roughly  equal  amounts.  Based  on  Fante’s  numerical  results,  about  half  the  turbulence  area  Ajis  beam 
spread  and  half  is  image  dancing. 


A  Possible  Scaling  Law  Parameter 

In  all  scenarios  in  which  the  laser  platform  velocity  is  constant  and  the  target  trajectory  is  a  straight  line, 
the  minimum  transverse  flow  speed  of  the  air  past  the  infinitely  extended  laser  beam  line  is  independent  of 
both  the  speed  and  the  location  of  the  target.  For  slewed  laser  beams,  the  proof  is  the  following:  In  the  rest 
frame  of  the  laser  platform,  the  orientation  of  the  plane  defined  by  the  laser  beam  and  the  missile  trajectory  is 
constant  and,  thus,  so  is  the  direction  of  its  normal.  Since  the  component  of  the  transverse  wind  velocity  in  the 
plane  of  the  laser  beam  and  the  missile  trajectory  is  always  cancelled  by  slewing  at  some  point  on  the  infinitely 
extended  laser  beam  line,  the  constant  normal  component  is  the  minimum  transverse  flow  speed.  Although 
deflection  of  the  laser  beam  out  of  the  above  plane  (by  thermal  blooming)  destroys  this  exact  invariance  (if  the 
angle  varies  with  the  location  of  the  missile  along  its  trajectory),  typical  deflections  are  less  than  1  m  in  1  km; 
so  the  minimum  transverse  flow  speed  is  still  nearly  a  constant  of  the  motion,  and  it  should  be  a  useful 
parameter  for  scaling  laws. 


Summary 

In  the  scenarios  investigated  here,  the  degradation  of  the  laser  beams  by  thermal  blooming  was  relatively 
modest  in  spite  of  the  stagnation  zone.  Only  for  the  CO2  laser  at  the  longest  range  and  the  shortest  offset  was 
the  beam  area  due  to  diffraction  and  thermal  blooming  larger  than  that  due  to  either  turbulence  or  jitter  alone. 
The  DF  laser  was  the  better  performer  of  the  two  lasers:  jitter  and  turbulence  dominated  thermal  blooming  at 
all  ranges  investigated. 
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6.  CIRCULAR  VS  ANNULAR  LASER  BEAMS 


The  initial  intensity  profile  of  a  high-energy  laser  beam  can  affect  the  on-target  focused-beam  intensity  by 
more  than  a  factor  of  2.  In  a  laboratory  experiment  that  combined  shape  change  with  optimal  cylindrical 
refocusing  for  each  shape,  Wallace  et  have  demonstrated  a  factor-of-5  increase  in  peak  on-target 

irradiance  by  replacing  a  Gaussian  beam  with  a  linear-ramp  intensity  profile  across  a  square  aperture,  one  of 
whose  sides  (as  well  as  the  incline  of  the  ramp)  was  aligned  with  the  wind.  Relative  to  a  Gaussian  beam  of 
equal  power  whose  e  diameter  matches  the  aperture  side  (or  diameter),  the  less  exotic  uniformly-illuminated 
square-  and  circular-aperture  laser  beams  seem  to  improve  the  on-target  irradiance  by  a  more  modest  factor  of 
2.31 

Understanding  the  effect  of  a  central  obscuration  is  especially  important  since  some  of  the  more  viable 
design  alternatives  employ  either  a  Cassegrain  telescope  (as  a  beam  expander  or  unstable  resonator)  or  a  coax¬ 
ial  pair  of  discharge  electrodes.  The  experimental  results  of  Phillips  et  suggest  that  obscuring  the  central 
10  to  15%  of  the  area  of  a  uniformly-illuminated  circular  aperture  (maintaining  the  power  level)  would  only 
decrease  the  on-target  irradiance  by  20  to  30%  when  the  thermal  blooming  is  strong.  However,  Gebhardt’s  3^ 
recent  calculations  suggest  that  decreases  of  up  to  at  least  67%  are  possible  in  specific  Army  scenarios.  We 
have  repeated  Gebhardt’s  calculations  on  our  FOUR-D  code  and  have  verified  his  conclusions. 


Table  7.  Circular  vs  annular  laser  beams;  target  plane  areas. 


Normalized  Beam  areas^ 


Laser 

Absorption 

No. 

Fresnel 

No. 

Slewing 

No. 

N 

Distortion 

No. 

FOUR-D 

code 

Gebhardt 

Annulus 

Circle 

Annulus 

Circle 

CO, 

0.40 

32.5 

74 

464 

3.0 

=  1.0 

4.2 

1.4 

CO, 

0.74 

17.6 

40 

856 

29.9 

16.2 

60.2 

34.2 

DF 

0.25 

49.2 

40 

596 

5.0 

2.3 

^Smallest 

area  containing  (1 

-  e'^l)  of  the 

target  plane 

power.  The  average 

intensity  is 

computed  for  this  area. 

Figure  22,  Initial  steady-state  heated^air  density  profiles  for  uniformly  illuminated  (a)  annular  (diameter  ratio  =  0.365)  and  (b)  circular 
aperture  cw  or  high-overlap-number  multipulse  laser  beams. 
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The  uniformly-illuminated  circular  and  annular  aperture  laser  beams  had  equal  power  levels  and  equal 
outer  diameters.  The  ratio  of  the  inner  to  the  outer  diameter  of  the  annulus  was  0.365.  The  performance  of 
both  laser  beams  was  evaluated  at  two  ranges  for  an  incoming-target  scenario  with  a  constant  overlap  number 
of  35  at  the  laser  aperture.  The  remaining  parameters  are  summarized  in  Table  7  in  dimensionless  form. 
Steady-state  thermal  blooming  was  assumed,  and  t  ^-blooming  was  ignored  because  the  saturation  times  ^  of 
63  MS  (2  km)  and  71  ms  (4  km)  are  much  longer  than  the  10-ms  pulse  duration. 


Annular  aperture 
(diameter  ratio  =  0.365) 


Circular  aperture 


Figure  23.  Isointensity  contour  plots  comparing  the  evolution  of  a  multipuise  CO2  laser  beam  with  uniformiy-illuminated  circular  aperture 
to  one  with  a  uniformiy-iliuminated  annular  aperture.  Both  iaser  beams  have  the  same  outer  diameter.  The  ratio  of  the  inner  to  the  outer  an¬ 
nular  diameter  is  0365.  The  dimensionless  parameters  of  this  case  are  =  0.74,  Nf  =  17.6,  =  40,  =  856,  and  Nq  =  35. 
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Our  results  are  summarized  in  Table  7,  along  with  Gebhardt’s  on-target  beam  areas.  Although  our 
on-target  beam  areas  are  significantly  smaller  (30%  smaller  for  N f=  32.5  and  50%  smaller  for  Nf=  17.6)  and 
our  intensities  are  higher  than  Gebhardt’s,  the  ratio  of  the  on-target  area  of  the  obscured  to  that  of  the  unob¬ 
scured  laser  beam  agree  within  10%  for  both  CO2  laser  calculations.  The  DF  laser  results  in  Table  7  are  for  75% 
of  the  power  and  for  the  same  scenario  and  range  as  the  longer  range  CO2  laser  calculations.  The  ratio  of  the 
target-plane  area  of  the  annular  beam  to  that  of  the  circular  beam  is  2.2  for  the  DF  laser,  compared  to  1.9  for 
the  CO2  laser. 

For  cw  or  high-overlap  multipulse  laser  beams,  the  air-density  profiles  due  to  laser  heating  at  the  aper¬ 
ture  are  shown  in  Fig.  22.  Along  any  line  parallel  to  the  wind,  the  density  profile  created  by  the  uniformly- 
illuminated  circular  aperture  is  a  linear  ramp.  Along  any  line  perpendicular  to  the  wind,  it  is  a  continuous 
curve  with  a  positive  curvature.  The  linear  ramp  deflects  the  laser  beam  into  the  wind  while  the  positive  curva¬ 
ture  defocuses  it  in  the  orthogonal  direction.  The  density  profile  for  the  annular  beam  is  simply  the  difference 
between  the  density  profiles  for  two  circular  apertures.  Along  a  line  perpendicular  to  the  wind,  the  density 
profile  in  the  shadow  of  the  depression  on  the  down-wind  side  of  the  annulus  has  negative  curvature.  This 
negative  curvature  causes  a  local  focusing  of  the  laser  beam,  as  shown  in  Fig.  23.  For  the  distortion  and  slew¬ 
ing  numbers  of  these  calculations,  the  focal  plane  is  clearly  far  in  front  of  the  target  plane;  thus,  subsequent 
diffraction  and  enhanced  thermal  blooming  due  to  the  higher  intensity  together  degrade  the  performance  of 
the  annular  beam  relative  to  that  of  the  circular  beam. 


7.  HIGH-ENERGY  LASER  BEAM 
PROPAGATION  THROUGH  A  SMOKE  SCREEN 


Gebhardt  has  provided  the  parameters  of  a  typical  white-phosphorous  smoke  cloud  and  of  an  engage¬ 
ment  in  which  such  a  cloud  might  occur.  The  geometry  of  this  engagement  is  described  in  Fig.  24  (unless  indi¬ 
cated  otherwise  the  smoke  source  is  located  at  Sj).  Since  the  vertical  dimension  of  the  smoke  cloud  and  the  dis¬ 
tance  between  the  laser  beam  and  the  smoke  source  are  large  compared  with  the  laser-beam  diameter,  the 
cloud  was  treated  as  vertically  infinite  and  the  variation  of  its  thickness  across  the  laser  beam  was  ignored.  The 
laser  beam,  target,  and  smoke  cloud  source  were  all  coplanar  and  the  target  and  wind  velocities  were  in  this 
plane.  The  smoke-cloud,  engagement,  and  laser  beam  parameters  are  tabulated  in  Table  8  as  dimensionless 
numbers.  Equal-diameter,  equal-power  CO2  and  DF  laser  beams  were  used;  they  had  a  uniformly-illuminated 
circular  aperture  beam  shape  and  were  focused  on  the  target. 

The  results  are  summarized  in  Table  9.  For  the  CO2  laser  the  space-averaged  on-target  intensity  is  re¬ 
duced  by  factors  of  4.4  and  21,  due  to  the  smoke  cloud  for  day  and  night  conditions,  respectively.  Of  this 
reduction,  only  a  factor  of  2  is  due  to  thermal  blooming  in  the  smoke  cloud;  the  remainder  is  due  to  unavoid¬ 
able  linear  absorption. 

The  effect  of  the  smoke  cloud  is  more  drastic  for  the  DF  laser,  due  to  the  smaller  atmospheric  absorption 
coefficient  and  smaller  diffraction-limited  focal-spot  size  at  this  wavelength.  The  on-target  space-averaged 
intensity  is  reduced  by  factors  of  26  and  320  for  day  and  night  conditions,  respectively.  A  factor  of  5  in  each  of 
these  reductions  is  due  to  thermal  blooming,  so  that  linear  absorption  and  thermal  blooming  have  comparable 
effects  during  the  day,  while  linear  absorption  dominates  at  night. 
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Compared  with  the  CO2  laser  in  the  same  scenario,  the  DF  laser  delivers  less  energy  to  the  target,  but 
delivers  it  to  a  smaller  area;  with  the  daytime  smoke  cloud  the  average  irradiance  of  the  DF  laser  is  20%  higher 
than  that  of  the  CO2  laser,  but  with  the  nighttime  smoke  cloud  it  is  30%  lower. 

The  DF  on-target  irradiance  for  daytime  conditions  with  and  without  the  smoke  cloud  is  shown  in 
Fig.  25.  Without  the  smoke  cloud  [Fig.  25(a)],  the  focal  spot  is  only  slightly  elongated  perpendicular  to  the 
wind  with  a  major-to-minor  diameter  ratio  (at  the  30%-contour  level)  of  1.4  and  is  deflected  less  than  one 
minor  diameter  downwind.  With  the  smoke  cloud  [Fig.  25(b)],  this  major-to-minor  diameter  ratio  grows  to 
6.7,  and  the  focal  spot  is  deflected  about  four  minor  diameters  downwind. 


Figure  24.  Scenario  for  smoke  cloud  calculations.  The  laser,  target,  and  smoke  source  are  at  the  points  L,  T,  and  Sj  (or  S2  ),  respectively. 


Table  8>  Smoke  cloud  propagation  conditions. 


Parameter 

DF 

Wavelengdi 

S 

8 

Fresnel  No. 

93.5 

33.5 

Slewing  No. 

101.0 

101.0 

Distortion  No.  (w/o  cloud) 

245.0 

263.0 

Absorption  No.  (w/o  cloud) 

0.170 

0.509 

Smoke  cloud 

Width  (a  ) 

Day 

9.19  X  10"^  R 

9.19  X  10"^  R 

y 

Ni^t 

3.89  X  10“^  R 

3.89  X  10"^  R 

Absorption  No. 

Day 

1.26 

1.01 

Night 

2.90 

2.33 

Scattering  No.^ 

Day 

0.473 

0.018 

Night 

1.09 

0.043 

^Scattering  coefficient  times 

cloud  thickness  (2<7^)* 

Table  9.  Beam  degradation  due  to  smoke  cloud. 

Laser  and 

Cloud 

Ratio  of 

Normalized 

smoke  cloud 

Inverse  of 

transmission 

average 

half-power  area 

condition 

area  ratio 

coefficient 

intensities 

^\I2 

CO^  day 

0.63 

0.357 

0.0321 

7.9 

COj  day® 

0.50 

0.357 

0.0255 

10.0*> 

COj  night 

0.51 

0.0935 

0.0068 

9.8 

CO^  w/o  smoke 

-  1.0 

1.0 

0.142 

5.0 

DF  day 

0.22 

0.177 

0.0389 

4.6 

DF  day® 

0.13 

0.177 

0.023 

8.0 

DF  ni^t 

0.17 

0.0185 

0.0032 

5.8 

DF  w/o  smoke 

=  1.0 

1.0 

=  1.0 

=  1.0 

^Cloud  positioned  at  S2  in  Fig.  19. 

^The  normalized  half-power  area  is  61.2  due  to  a  broad  low-intensity  wing  of  the  on-target  irradiance;  the  tabulated  number,  which 
corresponds  to  30%  of  the  power,  give  a  better  indication  of  the  performance  of  the  CO2  laser  in  this  scenario. 


Figure  25.  On-target  isointensity  contours  for  the  DF  laser  (a)  without  a  smoke  cloud  and  (b)  with  a  smoke  cloud  under  daytime  conditions. 
The  reference  line  indicates  the  relative  scales  of  the  two  plots. 
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8.  COMPARISON  OF  PROPAGATION  CONDITIONS 
AT  HELSTF  AND  NORTH  OBSCURA  PEAK 


Calculations  were  performed  to  evaluate  the  relative  merits  of  two  potential  laser  propagation  test  sites  at 
the  White  Sands  Missile  Range.  The  meteorological,  scenario,  and  laser  beam  parameters  for  the  calculations 
were  provided  by  Gail  Bingham,  who  also  evaluated  Atmospheric  Sciences  Laboratory  turbulence  data  and 
found  little  difference  between  the  two  sites.  We  have,  therefore,  neglected  turbulence  in  our  calculations. 
Four  scenarios  were  selected  for  the  site  comparison: 

•  A  pure  ground-to-ground  test  (case  1)  with  the  laser  beam  perpendicular  to  the  prevailing  wind; 

•  A  pseudo-ground-to-ground  test  (case  IB)  with  the  target  mounted  on  a  sled  moving  at  60  mph  into 
the  prevailing  wind; 

•  An  incoming-target  test  (case  2)  with  a  constant-altitude  target  at  a  range  equal  to  5  times  the 
horizontal  offset  of  the  linear  trajectory.  (The  trajectory  was  perpendicular  to  the  horizontal  component  of  the 
wind;  the  altitude  above  the  laser  was  0.305  times  the  offset  at  HELSPO  and  0  at  North  Obscura  Peak.); 

•  A  crossing-target  test  (case  3)  with  the  target  flying  into  the  horizontal  component  of  the  wind  at  a 
range  equal  to  the  distance  of  the  closest  approach. 

The  dimensionless  parameters  for  these  scenarios  are  given  in  Table  10.  For  the  North  Obscura  Peak  site, 
the  wind  field  is  a  function  of  the  distance  from  the  cliff;  we  approximated  this  dependence  by  =  -3.5  X  (1  - 
z)  m /s,  =  0.707  X  (3.5  +  6.5z)  m  js,  and  =  v^,  where  z  is  the  distance  from  the  cliff  in  km  and  -y,  is  vertical 

(the  zenith). 

The  results  of  the  calculations  are  summarized  in  Table  11  and  the  corresponding  on-target  intensity 
profiles  (including  jitter)  are  displayed  in  Figs.  26  and  27.  For  the  incoming-  and  crossing-target  scenarios 
there  is  little  difference  between  the  on-target  intensity  profiles  at  the  two  sites  because  jitter  is  the  dominant 
cause  of  beam  degradation.  If  the  jitter  were  ~3  /lirad/axis,  then  the  North  Obscura  Peak  Site  would  have  a 
slight  advantage  for  the  incoming-target  scenario.  The  maximum  potential  advantage  is  given  by  the  inverse  of 
the  ratio  of  the  jitterless  beam  areas  in  Table  11,  corrected  for  the  different  absorption  numbers:  1.6  for  the 
DF  laser  and  3.4  for  the  CO2  laser. 

For  ground-to-ground  tests  in  the  rainy  season,  the  propagation  conditions  at  North  Obscura  Peak  are 
significantly  better  than  at  the  HELSTF  site.  However,  addition  of  a  60-mph  target  sled  on  a  railroad  track  to 
the  HELSTF  facility  would  make  the  propagation  conditions  using  it  better  than  propagation  conditions  for 
stationary  target  ground-to-ground  testing  at  North  Obscura  Peak,  even  in  the  rainy  season.  Propagation  con¬ 
ditions  that  include  more  than  3  ^rad  of  jitter  do  not  dictate  a  choice  between  the  HELSTF  and  North 
Obscura  Peak  sites. 


Table  10.  Propagation  parameters  for  site  evaluation. 


Laser  Scenario 

HELSTF 

N.  Obscura  Peak 

Absorption 

No. 

Slewing 

No. 

Distortion 

No. 

Overlap 

No. 

Absorption 

No. 

Slewing 

No. 

Distortion 

No. 

Overlap 

No. 

DF®  1 

0.074 

0 

224 

_ 

0.055 

0 

74 

1  (dry) 

0.039 

0 

118 

- 

0.037 

0 

50 

- 

1  B 

0.074 

12.2 

224 

- 

- 

- 

- 

- 

2 

0.074 

30.1 

224 

- 

0.055 

13.7 

74 

- 

3 

0.074 

142.0 

224 

- 

0.055 

63.0 

74 

- 

C02^  1 

0.2 

0 

217 

36.5 

0.12 

0 

58 

16.2 

1  (dry) 

0.1 

0 

108 

14.3 

0.1 

0 

48 

16.2 

1  B 

0.2 

12.2 

217 

36.5 

- 

- 

- 

- 

2 

0.2 

30.1 

217 

35.6 

0.12 

13.7 

58 

16.3 

3 

0.2 

142.0 

217 

36.5 

0.12 

63.0 

58 

22.9 

^For  all  the  scenarios  a  fixed  range  and  beam  diameter  were  used.  The  Fresnel  numbers  were  265  for  the  DF  laser  and  95  for  the 
CO2  laser,  A  low  frequency  pointing-and-tracking  jitter  of  10  jurad  (rms/axis)  was  used;  turbulence  was  neglected.  A  uniformly- 
illuminated  annulus  with  an  inner-  to  outer-diameter  ratio  of  0.371  was  used. 
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Table  11.  Results  of  site  evaluation  calculations. 


Laser 

Scenario^ 

Normalized 
With  jitter 

beam  areas'’ 

Without  jitter 

Normalized 
average  intensity^ 

HELSTF 

N.  Obscura 
peak 

HELSTF 

N.  Obscura 
peak 

HELSTF 

N.  Obscura 
peak 

DF 

1 

<23.0 

^13.0 

- 

<2.7 

>0.037 

^0.092 

1  (dry) 

5.2 

3.5 

3.6 

0.49 

0.24 

0.23 

1  B 

2.0 

0.18 

- 

0.49 

- 

2 

1.3 

1.04 

0.052 

0.033 

0.81 

1.02 

3 

-  1.0 

0.85 

0.023 

0.028 

=  1.0 

1.3 

^^2 

1 

<42.0 

7.9 

2.9 

>0.014 

0.10 

1  (dry) 

7.0 

6.3 

6.0 

1.7 

0.11 

0.13 

1  B 

3.6 

- 

1.2 

- 

0.23 

- 

2 

2.0 

1.5 

0.57 

0.18 

0.37 

0.40 

3 

1.6 

1.4 

0.15 

0.15 

0.31 

0.38 

^All  values  are  for  the  rainy  season  unless  otherwise  indicated. 
^Area  within  the  31%  intensity  contour. 

^Intensities  include  the  effect  of  jitter  and  are  time-averaged  values. 
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Figure  26.  On-target  isointensity  contours  for  CO  2  laser  beams  for  White  Sands  Missile  Range  laser  test  site  selection  calculations.  These 
contour  plots  are  all  drawn  to  the  same  scale. 
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Figure  27.  On-target  isointensity  contours  for  DF  laser  beams  for  White  Sands  Missile  Range  laser  test  site  selection.  Contour  plots  are  all 
drawn  to  the  same  scale  used  in  Fig.  26. 
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